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Abstract. The nucleus-nucleus potential and energy dissipation in fusion reactions are obtained from micro- 
scopic mean-field dynamics. The deduced potentials nicely reproduce the one extracted from experimental data. 
Energy dissipation shows a universal behaviour between different reactions. Also, the dispersion of mass distri- 
bution in transfer reaction is investigated in a stochastic mean-field dynamics. By including initial fluctuations 
in collective space, the description of the dispersion is much improved compared to that of mean field only. The 
result is consistent with the macroscopic phenomenological analysis of the experimental data. 



1 Introduction 

The interplay between nuclear structure and dynamical ef- 
fects is crucial to properly describing nuclear reactions at 
energies close to the Coulomb barrier. Therefore, the the- 
ories describing such nuclear reactions need a unified de- 
scription for both nuclear structure and dynamics. More- 
over, recent developments on nuclear facilities introduce 
much interest on the properties of nuclei far from the sta- 
bility. The time-dependent Hartree-Fock (TDHF) theory 
H1I2I3I41 based on the Skyrme energy density functional 
(EDF) provides a rather unique tool for describing nuclei 
over the whole nuclear chart. The TDHF theory solves the 
time evolution of single-particle wave functions according 
to 

ifi-p = [h[p], P ], (1) 
ot 

where h[p] denotes the self-consistent mean-field Hamil- 
tonian from the Skyrme EDF, denoted by &[p], obtained 
from h[p] = 6&[p]/5p with the one-body density p. This 
model automatically includes important dynamical effects 
such as vibrations of nuclei, neck formations, and nucleon 
transfer during reactions. Since recent computational de- 
velopments now enable us to include all the terms of the 
Skyrme EDF used in static Hartree-Fock calculations in 
the three-dimensional coordinate space |5 6 7 8|, the de- 
scription of nuclear reactions using TDHF should be revis- 
ited. 

In this contribution, as an illustration of applications of 
the TDHF theory to nuclear reactions, we discuss the prop- 
erties of nucleus-nucleus potential and energy dissipation 
extracted from the TDHF model [9 101 . Moreover, we in- 
vestigate fluctuations of one-body observables, especially, 
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the dispersion of mass distributions in transfer reactions 
using a stochastic mean-field model 111 11121 . 

2 Nucleus-nucleus potential and energy 
dissipation from mean-field dynamics 

Nucleus-nucleus potential and energy dissipation are ex- 
tracted as follows 1191101 : (i) The TDHF equation for head- 
on collision is solved to obtain the time evolution of the 
total density of colliding nuclei, (ii) From the total den- 
sity, we compute at each time the relative distance R, as- 
sociated momentum P, and reduced mass /j of colliding 
nuclei, (iii) We assume that the time evolutions of R and 
P obey a classical equation of motion including a friction 
term which depends linearly on the velocity R: 

dR P dP dV 
dt /j. dt dR 

where V(R) and y(R) denote the nucleus-nucleus poten- 
tial and friction coefficient, respectively. The friction coef- 
ficient y(R) describes the effect of energy dissipation from 
the macroscopic collective degrees of freedom to the mi- 
croscopic ones. For the TDHF calculations presented in 
this contribution, the three-dimensional TDHF code devel- 
oped by P. Bonche and coworkers with the SLy4d Skyrme 
effective force |5| is used. The mesh sizes in space and in 
time are 0.8 fm and 0.45 fm/c, respectively. For more de- 
tails, see Refs. |f9TT0Tl . 

Dynamical effect on potentials deduced from TDHF 
trajectories at center-of-mass energies close to the Coulomb 
barrier is seen in all reactions considered here. Figure Q] 
shows the difference between the barrier height deduced 
from TDHF trajectories (Vb) and the barrier height ex- 
tracted from experimental data (V^ xp ) H 1 31 1 411 as a function 
of Vb- The solid line corresponds to the barrier height ex- 
tracted using high-energy TDHF trajectories (£ c . m . » Vb), 
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Fig. 1. Barrier height V B deduced from our method minus ex- 
perimental barrier height V| xp as a function of extracted barrier 
height for the reactions indicated in the figure. The values V B 
are deduced from high energy (solid line) and from low energy 
(dashed line) TDHF trajectories, respectively. 



whereas the dashed line is the result for low-energy TDHF 
trajectories (£ CJ1 ~ Vb)- The former identifies with the 
barrier height of the frozen density apploximation [15|. 
Dynamical reduction of the barrier height from high-energy 
TDHF to low-energy TDHF is clearly seen for all reac- 
tions. Moreover, because of this reduction, the value of the 
barrier height at low energy approaches the experimental 
data. This underlines the importance of dynamical effects 
close to the Coulomb barrier and shows the precision of 
our method. 

Our method is also able to provide information on en- 
ergy dissipation through the friction coefficient y, which is 
shown in Fig. [2] In this figure, we present reduced friction 
coefficients f3(R) = y{R)lp as a function of R scaled by 
the Coulomb barrier radius Rb for different reactions. Fig- 
ure|2]clearly shows that the order of magnitude of fi(R) and 
the radial dependence are almost independent on the size 
and asymmetry of the system. We also compare our results 
with that of a microscopic model based on small amplitude 
response by Adamian et al. (16) by the solid circles. The 
radial dependence and the magnitude of the friction coeffi- 
cient are very similar to those extracted from our method. 

From here, we conclude that the use of the macroscopic 
equation (0 is valid and mean-field dynamics gives good 
descriptions for the nucleus-nucleus potential and energy 
dissipation. 

3 Mean-field fluctuations 

It is well known that the mean-field model gives good de- 
scriptions for average evolution of one-body observables 
in low-energy nuclear reactions. However, it completely 
fails in the description of the dynamics of fluctuations of 
one-body observables. One of the shortcomings may be the 
failure of the mean-field description of the dispersion of 
mass distributions of final fragments in deep inelastic col- 
lisions. It has been recognized for a long time that TDHF 
calculations severely underestimate the dispersion of mass 
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Fig. 2. Extracted reduced friction j}(R) = y(R)/fi as a function of 
R scaled by the Coulomb barrier radius R B for different reactions. 
A microscopic friction coefficient by Adamian et al. is shown by 
the solid circles for comparison. 



distributions of experiments II17I18I . although TDHF cal- 
culations well reproduce the mean value of fragment mass. 

During the past decades, much effort has been devoted 
to overcoming this difficulty and to developing transport 
theories that are able to describe not only mean values but 
also fluctuations (for a review, see Refs [19 20]). Among 
them, the variational principle by Balian and Veneroni (B V) 
appears as one of the most promising methods B21I22I23I . 
However, even nowadays it remains difficult to apply 1 24 1 . 
More than 30 years after the first application of the TDHF 
theory, the absence of a practical solution to include fluctu- 
ations beyond mean field in a fully microscopic framework 
strongly restricts applications of mean-field-based theories. 

In order to overcome this difficulty, recently we pro- 
posed a stochastic mean-field (SMF) approach, which is a 
stochastic extension of the mean-field model for low en- 
ergy nuclear dynamics so as to include zero-point fluctua- 
tions of the initial state [11 25 1. The initial density fluctua- 
tions are simulated by representing the initial state in terms 
of a suitable ensemble of initial single-particle density ma- 
trices, which is similar to the idea in Refs. [26 27|. In fact, 
this idea can be regarded as the beginning of constructing 
time-dependent version of configuration mixing calcula- 
tions. In this manner, the description with single Slater de- 
terminant is replaced by a superposition of multi Slater de- 
terminants. A member of the ensemble, indicated by event 
label A, can be expressed as 

pV, r', = 2 <PW, ^)p-/cTT)<Zy T (r', t; A), (3) 

ijCTT 

where summations i and j run over a complete set of single- 
particle wave functions <£ 1(TT (r, f; A), and <x and t denote 
spin and isospin quantum numbers. According to the de- 
scription of the SMF approach 11251 . the elements of den- 
sity matrices p^.{crf) are assumed to be time-independent 

random Gaussian numbers with mean value pf(cn-) = 5ijn^ T 
and with the variance of the fluctuating part SpyXor) spec- 
ified by 

6pf j {o-T)8p A j , i ,{(r>T') 
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Here, n°~ T denotes the average single-particle occupation 
factor. At zero temperature occupation factors are and 
1, and at finite temperature they are determined by the 
Fermi-Dirac distribution. The great advantage of the SMF 
approach is that each Slater determinant evolves indepen- 
dently from each other following the time evolution of its 
single-particle wave functions in its self-consistent mean- 
field Hamiltonian, denoted by h(p A ), according to 

iH-& i(TT (r, t; A) = h(p A )& ilTT (r, t; A). (5) 
ot 

In the following applications, we focus on the head-on 
40 Ca+ 40 Ca collision around the Coulomb barrier energy. 



3.2 Mass dispersion in transfer reactions 

Another application has been made to the dispersion of the 
fragment mass distribution to improve the mean-field de- 
scription [12]. To do so, we investigate transfer reactions 
near the Coulomb barrier, where nucleon exchange will 
occur during reaction, and estimate the dispersion of the 
fragment mass distribution. 

In the SMF approach, time evolution of the mass num- 
ber of the projectile-like fragment A p is also described by 
a Langevin equation |29|, 



-A A =v(A P , f) +£(*), 



(8) 



where v(A A , f) denotes the drift term for nucleon transfer. 
The Gaussian white noise random force £j(0 is determined 
with zero mean value and with a correlation function, 



3.1 Fusion reactions 



e A (t)e A (t') = 2S(t-t')D A 



(9) 



First, we apply the SMF approach to fusion reactions ifTTl . 
To discuss the fluctuation of collective variables, we map 
the SMF time evolution to a one-dimensional macroscopic 
Langevin equation, which is similar to Eq. (O except an 
additional Gaussian random force g p (t): 



d_ p A 
dt 



■ y(R A )R A + #(f), 



(6) 



Ignoring non-Markovian effects, the random force g p (f) 
with zero mean value reduces to white noise specified by a 
correlation function, 



= 2S(t - t')D PP (R). 



(7) 



Here D PP (R) denotes the momentum diffusion coefficient. 
We note that the expression of the diffusion coefficient has 
the same form as that obtained from the phenomenologi- 
cal nucleon exchange model [28]. As an example, the dif- 
fusion coefficient for the head-on 40 Ca+ 40 Ca collision at 
£ c m = 100 MeV is shown in Fig. [3] 
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Fig. 3. Diffusion coefficient as a function of the relative distance 
for the head-on 40 Ca+ 40 Ca collision at E c m = 100 MeV. 



where D AA is the diffusion coefficient associated with nu- 
cleon exchange. The variance cr 2 AA of fragment mass dis- 
tribution is determined by small fluctuations of the mass 

number 5A A p through o~ AA (t) = 6A p 6A p . According to the 
Langevin equation, neglecting contributions from the drift 
term, the variance is related to the diffusion coefficient ac- 
cording to H29I301 



'AA 



(t) 



-2 f D A 
Jo 



(10) 



In the phenomenological nucleon exchange model, the 
relation cr 2 AA (f) - N exc (t) was obtained, where N exc (t) de- 
notes the accumulated total number of exchanged nucleons 
until time f, and was extensively used to analyze the exper- 
imental data of mass dispersion ||3"D . In the following, to 
check whether the SMF approach satisfies the above rela- 
tion, we estimate the both quantities by the SMF approach. 

In Fig. |4] the variances of the fragment mass distri- 
butions deduced from the SMF approach for the head-on 
40 Ca+ 40 Ca reaction at three center-of-mass energies are 
shown by lines. The number of exchanged nucleons is in- 
serted in Fig. [4] by the solid circles, solid triangles, and 
solid squares from high to low energies. The mass vari- 
ance estimated from the SMF approach is consistent with 
this relation. We also estimate the variance of fragment 
mass distribution using the standard TDHF approach. The 
asymptotic values of cr AA for the 40 Ca+ 40 Ca reaction are 
0.004, 0.008, and 0.008 from low to high energies, while 
the number of exchanged nucleons are 0.432, 1.441, and 
3.634. The TDHF results are much smaller than the num- 
ber of exchanged nucleons and are also much smaller than 
the results obtained from the SMF approach that are 0.730, 
1.718, 3.790. The failure of the TDHF theory on the de- 
scription of variances of the fragment mass distribution has 
been recognized for a long time as a major limitation of the 
mean-field theory. It appears that the SMF approach cures 
this shortcoming. As seen from Fig.|4] not only the asymp- 
totic value of cr 2 AA but also the entire time evolution is very 
close to the evolution of N exc (t). 
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Fig. 4. Time evolution of cr^ obtained from the SMF ap- 
proach for 40 Ca+ 40 Ca reaction at different center-of-mass ener- 
gies. Number of exchanged nucleon is superimposed by the solid 
circles, solid triangles, and solid squares from high to low ener- 
gies. 



4 Conclusion 

Mean-field dynamics and mean-field fluctuations using mi- 
croscopic time-dependent models are discussed in the con- 
text of low energy nuclear reactions. We have shown that 
the TDHF theory gives precise values of nucleus-nucleus 
potential and a universal behavior of energy dissipation. 
We have also shown that the SMF approach correctly de- 
scribes the mass dispersion of final fragments in transfer 
reaction at energies near the Coulomb barrier. This gives a 
practical solution to properly describe mean-field fluctua- 
tions on top of mean field. 

This work is supported in part by US DOE Grant DE-FG05- 
89ER40530. 
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